Self-induced oscillations in an optomechanical system 
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We have explored the nonhnear dynamics of an optomechanical system consisting of an illuminated 
Fabry-Perot cavity, one of whose end-mirrors is attached to a vibrating cantilever. Such a system 
can experience negative light-induced damping and enter a regime of self-induced oscillations. We 
present a systematic experimental and theoretical study of the ensuing attractor diagram describing 
the nonlinear dynamics, in an experimental setup where the oscillation amplitude becomes large, 
and the mirror motion is influenced by several optical modes. A theory has been developed that 
yields detailed quantitative agreement with experimental results. This includes the observation of 
a regime where two mechanical modes of the cantilever are excited simultaneously. 



Micro- and nanomechanical systems have become a 
focus of both theoretical and experimental research [IJ, 
with the goals ranging from ultrasensitive measurements 
to fundamental tests of quantum mechanics. One partic- 
ularly promising branch of this field deals with optome- 
chanical systems, where the interaction of light (stored 
inside an optical cavity) with macroscopic mechanical de- 
grees of freedom (such as the coordinate of a cantilever) 
is exploited. This can give rise to a variety of effects, 
including a modification of the mechanical spring con- 
stant m 131 m [51 bistability [TJ [5], optomechanical 
cooling, and parametric instability. A recent series of 
experiments has demonstrated impressive progress with 
respect to coohng [51El|ll[ini[IIl[ia[I3[Il], which may 
ultimately lead to the quantum ground state [HI [16] of 
mechanical motion in such devices. On the other hand, 
the opposite regime is of equal interest, where the me- 
chanical Q factor is enhanced since the mechanical de- 
gree of freedom extracts energy provided by the optical 
radiation. In that regime, a parametric instability arises 
which drives the system into a state of self-sustained os- 
cillations [l[T71[ni[Il[ini[lIl[21[23[21!25]. More- 
over, it is now known that the physics of that regime 
also applies to other systems as diverse as an LC cir- 
cuit driven by a radio-frequency source |26j or a current- 
biased superconducting single-electron transistor coupled 
to a nanobeam [27J [25] • Although the basic instability 
has been observed by now in a number of experiments 
[m Uni [m HI [Ml US], it was recently realized theo- 
retically f23l that the nonlinear dynamics of this system 
can become highly nontrivial, leading to an intricate at- 
tractor diagram. Here we report on an experiment that 
traces this diagram, and combine it with a detailed the- 
oretical analysis and systematic quantitative comparison 
of theory and experiment. In addition, this comparison 
has enabled us to observe an unexpected feature, the si- 
multaneous excitation of several mechanical modes of the 




Figure 1: (Color online) The experimental setup. The light 
inside the optical cavity is focused onto a cantilever, where 
it exerts a force. Both the average transmitted light and its 
sidebands at the cantilever frequency are recorded ("Trans- 
mission" and "Amplitude", respectively). 



cantilever, leading to coupled nonlinear dynamics. 

Experimental setup. - To study these questions, we em- 
ploy the optomechanical setup displayed in Fig. [T] The 
light of a 633 nm monomode HeNe-laser is coupled into 
a single mode fiber and passes through a Faraday isola- 
tor (35 dB suppression). The fiber end inside a vacuum 
chamber (at 5-10"^ mbar) was polished and coated with a 
reflecting gold layer of 30 nm (yielding a theoretical reflec- 
tivity of 70%) to form the first cavity mirror. The sample 
is a gold-coated AFM cantilever acting as a micromirror, 
with length 223 ^m, thickness 470 nm, width 22 /im, 
spring constant K = 0.01 N/m, and a gold layer of 42 nm 
evaporated on one side only. The cantilever's fundamen- 
tal mechanical mode has a frequency of cui = 2tt ■ 8.7 kHz 
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and a damping rate of Fi = 30 Hz. A simulation of the 
silicon-gold bilayer system gave a reflectivity of 91% for 
a wavelength of 633 nm. The divergent beam coming out 
of the fiber is sent through a microscope setup consisting 
of two identical lenses, yielding a Gaussian focus on the 
sample with a 1 /e^-diameter of 6 ^m. The cantilever has 
been mounted on an xyz piezo stepper positioner block 
[30] , such that it can be placed at the microscope's focal 
point, which was chosen near the end of the cantilever, 
at about 3/4 of its length. The finesse of the cavity de- 
fined by the sample and the fiber end was found to be 
F 4.5. The transmitted intensity is measured with a 
Si photodiode behind the cantilever, while sweeping the 
cantilever position through the optical resonance. 

Theoretical model. - The dynamics of the cantilever is 
described by the equation of motion of a damped oscil- 
lator, driven by light-induced forces: 



Self-induced oscillations. - Time-retarded forces in- 
duce an effective optomechanical damping rate whose 
sign changes when passing through the resonance [21 IH 
[TSl [TBI [m E3] . When the full damping rate becomes neg- 
ative, the system undergoes a Hopf bifurcation and the 
cantilever settles into stable oscillations [53] which (for 
the parameters of interest here) are sinusoidal to a very 
good approximation: x{t) = x + Acos{ijJit). The non- 
linear dynamics can then be characterized by solving for 
the amplitude A and offset x. From these, it will be pos- 
sible to obtain the experimentally observed evolution of 
the light intensity /(t). In steady state, the average force 
and power input must be zero, i.e. (i) = and [xx] = 0, 
where (. . .) denotes the time-average. Inserting Eq. ([l]), 
we obtain the power balance equation: 



-uj'i{x-x^)-Tix+{F''''' + F''°')lmi (1) 
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Here x{t) is the cantilever deflection observed at the laser 
spot. For now, we focus on the motion of the first me- 
chanical mode (with effective mass mi, frequency uji, and 
damping rate Fi), and we will return to the influence of 
higher-order modes below. The cavity is assumed to be 
in resonance with the laser when x = 0, while the me- 
chanical equilibrium position is given by a;o (sometimes 
referred to as the 'detuning'). In the experiment, it is 
controlled by the piezo positioning system. The radia- 
tion pressure force F''^'^ is proportional to the power / 
circulating inside the cavity, F'^^'^/rrii — V ■ I. The bolo- 
metric force F^°^ arises due to light being absorbed, thus 
heating the cantilever that then deforms as a bimorph. 
It is enhanced by a factor A over f and is retarded 
due to the finite time of thermal conductance r: 



The radiation pressure does not contribute, (xl) — 0, 
since the intensity follows the motion instantaneously 
(thus xl is an antisymmetric function of time). The ex- 
pression {x0) may be simplified by introducing the first 
harmonic of the light intensity, Ii — ]p dt cos{uJit)I{t) 
(with the period of mechanical motion T =^ 2tt/u!i): 



-Auji ■ Ii 
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This yields the first of the equations needed to find {x, A): 



F^°^/mi ^ AV 



dt' 



p-it't')/^I(^t')=AVe{t). (2) 



Here 9 is proportional to the change in temperature 
brought about by absorption of light. 

We restrict the following discussion to the case of an 
optical ring down time that is much smaller than the pe- 
riod of cantilever oscillations (like in the present setup). 
Thus, the light intensity reacts instantaneously to the 
cantilever motion, I{t) — I[x{t)]. Due to the small fi- 
nesse, we have to employ the Airy function dependence 
describing a series of overlapping Fabry-Perot resonances: 
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(3) 



Here /max is the peak circulating power, which is propor- 
tional to the input power. In contrast to the theory in 
|23j . where a large optical ringdown time led to intensity 
oscillations, here the nonlinear dynamics is induced by 
the time retardation of the bolometric force. 



Fi —A uiT 
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dip I[x + A COS ip] COS ip (6) 



Moreover, the force balance condition 



ujl{x-x,)^V-{{I)+K{e)) (7) 



enables us to obtain x = x{x{), A): 



{l+K)'P f^'' 
X — xq — — ^_ 2 — / dtp I[x + A cos tf\. (8) 
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As a special case, for A = 0, this contains the physics of 
static optomechanical bistability [7J [H] . 
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Figure 2: (Color online). Theoretically predicted attractor di- 
agram for the nonlinear cantilever dynamics, (a) Amplitude 
A of optomechanically induced oscillations vs. equilibrium 
position xo of the cantilever. The colorscale plot displays the 
net power fed into the cantilever from the radiation field, i.e. 
the r.h.s. of Eq. |6] Contour lines indicate possible oscilla- 
tion amplitudes (dynamical attractors) for different values of 
the damping constant Fi. Red dots are the results of di- 
rect numerical simulations, for the value of Fi realized in the 
present experimental setup (and a laser input power level of 
/ — 0.3/o). The inset displays the same plot, but as a func- 
tion of {A,x) instead of (A, xq), i.e. without the distortion 
due to the time-averaged radiation force that leads to x xq. 
(b,c) Whenever the cantilever motion (thick black line) passes 
through the optical resonances (green bars), the light inten- 
sity (red line) displays spikes, which lead to delayed increases 
in the radiation force (dashed blue line), which was plotted as 
{0 — 9) X 200. Plots (b) and (c) correspond to the respective 
positions indicated in the attractor diagram (a). 

Attractor diagram. - Using Eqs. (|6| and ([8|, one ob- 
tains solutions (x, A) which can be visualized in attrac- 
tor diagrams, like the one shown in Fig. [2] The color 
scale encodes the power input due to the light-induced 



forces (r.h.s. of Eq. [6]), as a function of Xq and A. The 
solution of Eq. |6]for various values of Ti/V then corre- 
sponds to contour lines of this function. Apart from the 
expected A/2-periodicity in the detuning xq, the main 
feature is the appearance of multiple solutions for A at 
a given xq ("dynamical multistability" ). Self-induced 
oscillations appear when x{t) briefly dips into the res- 
onance near its turning point, thereby gaining energy. 
Thus, to a first approximation, attractors appear along 
the diagonals A w x + nX/2. Near the special points 
X = A/4 -t- n A/2, the total power input vanishes due to 
symmetry: The cantilever extracts energy from the light 
field at one turning point but loses an equal amount of 
energy at the other turning point. 

The deviation between x and Xq, obtained from Eq. 
([s]), leads to a distortion of the diagram (compare inset 
of Fig. |2]). This effect grows with increasing input power, 
finally leading to multiple solutions for x{xq, A). 

Comparison of theory and experiment. - In the ex- 
periment, the detuning xq and the input power 1^ are 
varied, while the transmitted light intensity is measured. 
This can be compared to the time-averaged circulating 
power {I{t)) obtained from the theory. Moreover, as soon 
as the self-induced oscillations set in, I{t) is modulated 
at the cantilever frequency. A very helpful nontrivial 
feature of this system is the existence of a relation be- 
tween the motional amplitude A and the first harmonic 
of the light intensity, Ii (defined above). From Eq. 
we see that they are directly proportional, with the pro- 
portionality factor containing only known, fixed param- 
eters: A = -2(AP/wir)u;iT(l + (wir)2)-i7i. This rela- 
tion holds true only in steady state (on the attractors), 
but then it is valid even when the motion sweeps across 
several optical resonances. Experimentally, the first har- 
monic Ii is obtained by sending the photodetector signal 
through a narrow bandpass filter (100 Hz) centered at the 
eigenfrequency uji of the first mechanical mode. 

Theoretical and experimental curves for both the av- 
erage intensity ("transmission") and the amplitude are 
shown in Fig. [3] for different input powers. For this 
comparison, the following parameters have been used: 
F « 4.5 (from a fit at low input power), A = 3950 and 
ujiT = 39 (obtained independently). The overall con- 
version factor between experimentally measured input 
power and the force on the cantilever was found to be 
VIo = 0.0775m/s^, by adjusting for a good fit to the 
data at intermediate power (the same was done for the 
rescaling of theoretical and experimental transmission in- 
tensity). Here the maximum laser power Iq = 1.3 mW is 
estimated to yield 500 /iW circulating in the cavity on 
resonance. 

At the lowest power displayed in Fig. [3] self-oscillations 
have just set in, and the transmission curve shows a strik- 
ing asymmetry. At higher input powers, the multistabil- 
ity predicted by the attractor diagram. Fig. |2] leads to 
hysteresis effects that appear upon sweeping xq up and 
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take into account the second mode as well. The total 
displacement is x{t) = Xo + xi{t) + X2(t), and we have to 
employ a set of equations: 



-A/4 A/4 -A/4 A/4 

Displacement x^, Displacement 

Figure 3: (Color online) Experiment vs. theory. The trans- 
mitted light intensity (left) and the amplitude of self-induced 
cantilever oscillations (right), from a simulation of the theo- 
retical model (red full curves) and from the experiment (blue 
data points), at increasing input power levels (top to bottom). 
Theoretical "transmission" curves display the (rescaled) time- 
averaged circulating light intensity from the simulation. "Am- 
plitude" curves are obtained from the power in the intensity 
sidebands at the fundamental mechanical mode of the can- 
tilever (see main text). For clarity, the hysteresis observed 
upon sweeping xo up or down has been shown only in the mid- 
dle panel. The region of instability (shaded interval) grows 
with increasing input power. Simultaneous self-induced oscil- 
lations of the first two mechanical modes set in at the highest 
power displayed (in the two intervals indicated in the plot). 
The calculated amplitude of the second mode is shown as a 
dashed line. 



down. 

Beginning at /in = 0.57 /g, a second interval of self- 
oscillatory behaviour appears to the left of the resonance, 
growing stronger and wider with increasing laser power. 
This initially completely unexpected result may be ex- 
plained by invoking the influence of higher mechanical 
modes. These may be excited by the radiation as well, 
leading to coupled (multimode) nonlinear dynamics with 
richer features than discussed up to now. 

In order to describe the behaviour in that case, we now 



-ujfxi 



F}'°'[x{t)]/m, 



(9) 



where Xi denotes the coordinate of the i-th mechanical 
mode with frequency uji, mechanical damping rate F;, 
and effective mass (where wi/27r — 8.7kHz, LOijli^ — 
60 kHz, Fi = 30.0 Hz, Fa = 150 Hz). We neglected the ra- 
diation pressure force, as this is much smaller anyway for 
the parameters of this setup. The mechanical modes are 
now coupled indirectly by the bolometric force. For the 
present setup, this force changes sign when going to the 
second mode. Choosing F^°^vn\l F^°^m2 — —28.8 as an 
adjustable parameter, we have found the numerical sim- 
ulation of these coupled nonlinear equations for the first 
two modes to be in surprisingly good agreement with the 
experiment (Fig. [3|. We have to note that the relation 
between the measured "amplitude" , i.e. first harmonic of 
/(i) at frequency uJi, and the actual amplitude does 
not hold exactly if both modes are excited simultane- 
ously. 

At maximum laser power, there are actually two inter- 
vals with simultaneous excitation of both modes (indi- 
cated in Fig. [3]). Specifically, the onset of such a regime 
at ~ A/8 can be interpreted as follows: Taking into ac- 
count F^°^ IF\°^ < 0, we see that the second mode gains 
its energy from dipping into the resonance at x — A/2, 
while the first is still provided with energy due to the 
resonance at a; = 0. 

Numerical evidence shows that the steady-state mo- 
tion consists of sinusoidal oscillations in 0:1.2 at the re- 
spective cigenfrequencies, of nearly constant amplitudes 
and without phase locking (for the parameters explored 
here). Thus x{t) fa xo + J2i=i[^i + cos(ijJit + (f>i), 

where SAi{t)/Ai <C 1, and the (j>i are arbitrary phases. 
Higher input powers will lead to excitations of additional 
modes, and the system might go into a chaotic regime of 
motion. 

Conclusions. - We have analyzed the nonlinear dy- 
namics of an optomechanical system, by measuring and 
explaining its attractor diagram. The comparison of 
data and theoretical predictions have revealed the on- 
set of multi-mode dynamics at large optical power, with 
two mechanical modes of the cantilever participating in 
the radiation-driven self-sustained oscillations. These ef- 
fects could find applications in highly sensitive force or 
displacement detection [53]. In the future, it would be 
interesting to observe the attractor diagram in systems 
of a high optical finesse [51] [22] (with delayed radiation 
dynamics) , or the self-excitation of multiple mechanical 
modes of sub- wavelength mechanical resonators interact- 
ing with the radiation field inside a cavity [HJ US]- The 
whole field of quantum nonlinear dynamics in systems of 
this kind also remains to be explored. 
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